Comparison of synthesized and acquired high b-value diffusion-weighted MRI for detection of prostate cancer

Background High b-value diffusion-weighted images (DWI) are used for detection of clinically significant prostate cancer (csPCa). This study qualitatively and quantitatively compares synthesized DWI (sDWI) to acquired (aDWI) for detection of csPCa. Methods One hundred fifty-one consecutive patients who underwent prostate MRI and biopsy were included in the study. Axial DWI with b = 0, 500, 1000, and 2000 s/mm2 using a 3T clinical scanner using a 32-channel phased-array body coil were acquired. We retrospectively synthesized DWI for b = 2000 s/mm2 via extrapolation based on mono-exponential decay, using b = 0 and b = 500 s/mm2 (sDWI500) and b = 0, b = 500 s/mm2, and b = 1000 s/mm2 (sDWI1000). Differences in signal intensity between sDWI and aDWI were evaluated within different regions of interest (prostate alone, prostate plus 5 mm, 30 mm and 70 mm margin and full field of view). The maximum DWI value within each ROI was evaluated for prediction of csPCa. Classification accuracy was compared to Restriction Spectrum Imaging restriction score (RSIrs), a previously validated biomarker based on multi-exponential DWI. Discrimination of csPCa was evaluated via area under the receiver operating characteristic curve (AUC). Results Within the prostate, mean ± standard deviation of percent mean differences between sDWI and aDWI signal were -46 ± 35% for sDWI1000 and -67 ± 24% for sDWI500. AUC for aDWI, sDWI500, sDWI1000, and RSIrs within the prostate 0.62[95% confidence interval: 0.53, 0.71], 0.63[0.54, 0.72], 0.65[0.56, 0.73] and 0.78[0.71, 0.86], respectively. Conclusion sDWI is qualitatively comparable to aDWI within the prostate. However, hyperintense artifacts are introduced with sDWI in the surrounding pelvic tissue that interfere with quantitative cancer detection and might mask metastases. In the prostate, RSIrs yields superior quantitative csPCa detection than sDWI or aDWI.


Background
Diffusion-weighted imaging (DWI) is a critical component of multiparametric MRI for the detection and characterization of clinically significant prostate cancer (csPCa) [1].The degree of diffusion-weighting in DWI is indicated by the b-value, with higher b-values corresponding to images with less signal where water in tissues diffuses more rapidly [2].High b-values are used for their greater tumor conspicuity and detection of even small lesions [3].The Prostate Imaging -Reporting and Data System (PI-RADS v2.1) recommends the acquisition of high b-values (1400-2000s/mm 2 ) for lesion detection, without precisely defining an optimal value for csPCa [4].While clinically valuable, high b-values require more scan time and suffer from low signal-to-noise ratio (SNR) and increased susceptibility to artifacts due to microscopic motion or small fluctuations in local magnetic field.One common solution, permitted by PI-RADS, is to synthesize high b-value images by extrapolating signal from acquired low b-value images using a mono-exponential model [1,5].However, mono-exponential models do not adequately represent restricted diffusion in complex tissues [6,7], possibly calling into question the accuracy of synthesized images.
Studies have yielded contradicting results on whether synthesized b-values are clinically interchangeable with acquired DWI (aDWI) images.Liu et al. [23] compared various models, including the standard mono-exponential, for the detection of csPCa and concluded that nonlinear fitting with various b-values is superior to simpler models.In contrast, other studies reported better image quality for synthetic DWI (sDWI) with a similar tumor detection rate in comparison to acquired DWI [5,[24][25][26][27].
In this study, we qualitatively and quantitatively analyzed the differences between acquired and synthesized high b-value images for detection of csPCa.Further, we evaluated acquired and synthesized high b-value DWI for detection of csPCa at the patient level.For comparison, we also evaluated RSIrs, a quantitative biomarker based on multi-compartment DWI that is known to perform well for patient-level csPCa detection.

Patient cohort
This retrospective study was approved by the institutional review board at UC San Diego (IRB 805394).The research was performed in accordance with the Declaration of Helsinki, and all relevant regulations.A waiver of consent was approved by the institutional review board for this study as there was minimal risk of harm to patients.The retrospective dataset was described previously [22].Briefly, 440 consecutive men who underwent prostate MRI examination with a multi-b-value diffusion acquisition (compatible with Restriction Spectrum Imaging, RSI) between November 2017 and December 2020 were considered for inclusion.Patients were excluded if they had undergone prior treatment for prostate cancer or if there was no available biopsy result performed within 180 days of MRI acquisition.In total patients were included in the study.Patient characteristics are summarized in Table 1.MRI examinations were interpreted per routine clinical practice by ten board-certified (median of four years of experience) and subspecialty fellowship-trained radiologists.For all patients, suspicious lesions were contoured per PI-RADS v2.1 using MIM software (MIM Software, Inc; Cleveland, OH).For the present study, whole-gland prostate segmentation was performed using OnQ Prostate software (Cortechs Labs, San Diego, CA, USA).Clinically significant prostate cancer (csPCa) was defined as grade group ≥ 2. In patients who underwent prostatectomy, grade group was determined per final pathology report.Biopsy (typically systematic and targeted) and prostatectomy were performed according to clinical routine, and both were examined by board-certified pathologists.86 of the 151 patients were found to have csPCa, while 65 had only benign tissue or grade group 1 cancer (further details in Table 1).

MRI acquisition
All MRI acquisitions were performed on a 3T clinical GE scanner (Discovery MR750, GE Healthcare, Waukesha, WI, USA) using a 32-channel phased-array body coil surrounding the pelvis.Acquisition parameters are summarized in Table 2.A single axial DWI volume was acquired for each patient.T 2 -weighted reference images were acquired for all patients with field of view (FOV) identical to the DWI volume.RSI calculations were performed as described in prior studies [20][21][22].
Post-processing of the image data was performed using in-house software in MATLAB (version R2017a, MathWorks, Natick, MA, USA).DWI images were corrected for B 0 inhomogeneity distortions, gradient nonlinearity, and eddy currents [28][29][30].Multiple acquired DWI samples at specific b-values were averaged together and normalized by median signal intensity of urine in the bladder at b = 0 s/mm 2 .

Synthetic b-value computation
Synthetic high b-value DWI (sDWI) was calculated using the conventional, mono-exponential formula (see below) and using b-values up to 500 s/mm 2 (sDWI 500 ) or b-values up to 1000 s/mm 2 (sDWI 1000 ).

S(b)
is DWI signal for a given b-value, b. S 0 is the signal with no diffusion weighting.ADC is the apparent diffusion coefficient.sDWI was calculated for b = 2000s/ mm 2 to match the acquired high b-value DWI (aDWI).To explore the application of sDWI and aDWI for detection of significant cancer lesions outside of the prostate, sDWI and RSIrs were additionally calculated for one representative patient with csPCa and bone metastasis.

Data analysis
All data analysis was performed using in-house MAT-LAB scripts (version R2021a, MathWorks, Natick, MA, USA).Quantitative differences between sDWI and aDWI were estimated by a voxel-wise comparison of the images.Relative deviations were calculated for three different regions of interest (ROIs): prostate, prostate plus a margin of 5 mm, and the whole field of view (FOV) using the following formula: where S s is the synthetic signal intensity, S a the acquired signal and N the number of voxels in the considered images.Mean and standard deviation of ΔS over all patients are reported.A negative value indicates that the acquired signal intensity is higher than the synthesized signal intensity.Further, violin plots were generated for the 50 th , 95 th , and 98 th percentile of signal intensity within several ROIs: prostate; prostate plus margin (5 mm, 30 mm, or 70 mm); and the whole FOV.For the whole FOV, values higher than 3000 signal intensity units (SIU) were capped and set to 3000 SIU.Violin plots present the median value in combination with the kernel density distribution [31].
Lesion conspicuity was evaluated using the contrast-tonoise ratio (CNR) between lesion and surrounding prostate tissue.CNR is defined as the following: where μ is the mean signal of the ROI under consideration and σ the standard deviation.CNR was evaluated for all patients and patients diagnosed with csPCa.A higher CNR indicates a better tumor conspicuity [32]. N/A Significant differences between CNRs of different images and patient cohorts were tested using two sample t-test with a confidence level of 0.01.Prediction of whether csPCa was found on biopsy was also evaluated for aDWI, sDWI, and RSIrs.RSIrs is a quantitative cancer biomarker based on a multi-exponential DWI model (four compartments) and has been previously shown to be more accurate than conventional DWI [20][21][22].Computation of RSIrs for this dataset was performed previously and is described in detail in previous publications [20][21][22].Briefly, the coefficient for the slowest diffusion compartment (corresponding to intracellular restricted diffusion) was normalized by the median signal within the prostate on b = 0 s/mm 2 images.The maximum aDWI, sDWI, or RSIrs value within each considered ROI was used as the predictor variable [22].This is analogous to the maximum standard uptake value (SUV) in quantitative Positron Emission Tomography (PET) imaging.Receiver-operating characteristic (ROC) curves were calculated, and the area under the curve (AUC) reported for aDWI, sDWI, and RSIrs.The false positive rate at 90% sensitivity (FPR90) was also reported for each metric to illustrate performance at one threshold [21].AUC and FPR90 were compared using bootstrap (N = 10,000) 95% confidence intervals and p-values.

Results
Figure 1 shows the difference between acquired and synthesized b-values for a representative patient using the mean signal intensity within the prostate.Within the prostate, mean ± standard deviation of percent differences between sDWI and aDWI were -46 ± 35% for sDWI 1000 and -67 ± 24% for sDWI 500 .A negative error indicates sDWI had lower intensity values than aDWI.Figure 2 shows aDWI, sDWI and RSIrs for three representative patients with the same window and level.
Comparing sDWI 1000 to sDWI 500 , a difference of -41 ± 4% was estimated (see Table 3).sDWI 500 had overall larger errors than sDWI 1000 .Signal intensity of sDWI was lower than aDWI in the prostate and in the prostate plus 5 mm margin considering all voxels, as indicated by a negative median difference, see Table 3. Further, a lower median value considering only the 50 th , 95 th and 98 th percentiles was also observed.The 50 th percentile of aDWI is higher than sDWI for all considered ROIs.For the 95 th and 98 th percentiles, however, sDWI is larger for margins of ≥ 30 mm beyond the prostate.The standard deviation of sDWI is larger than the standard deviation of aDWI for all considered percentiles in all ROIs.A comparison of the 50 th , 95 th and 98 th percentiles for five ROIs with varying margins around the prostate is shown in Fig. 3.
Figure 5 compares sDWI, aDWI and RSIrs for detection of significant cancer lesions outside of the prostate.For detection of csPCa, the AUCs for sDWI and aDWI were similar in both the prostate and prostate plus 5 mm (Fig. 4 and Table 4).Classification accuracy decreased significantly for sDWI when considering the whole FOV (AUC = 0.45 [0.36, 0.54] for sDWI 1000 and 0.47 [0.38, 0.56] for sDWI 500 ).RSIrs was superior to sDWI and aDWI for all ROIs (p < 0.01).The AUC of RSIrs was 0.77 [0.69, 0.84] within prostate plus 5 mm and decreased to 0.70 [0.61, 0.78] for the whole FOV.FPR90 was similar for aDWI and sDWI in all ROIs.Mean FPR90 was significantly lower for RSIrs than for either aDWI or sDWI, indicating fewer false positives (p < 0.05).

Discussion
We found that synthesized DWI images can be qualitatively similar to acquired DWI within the prostate even though sDWI is quantitatively an inaccurate representation of aDWI.Moreover, sDWI introduces unacceptable artifacts and inaccuracies in surrounding pelvic tissues.sDWI 500 calculated using only b-values up to 500 s/mm 2 was inferior to sDWI 1000 , demonstrating that the accuracy of sDWI is improved by including higher b-values.Acquiring additional b-values ≤ 1000 s/mm 2 would further increase the accuracy of the mono-exponential fit and the calculation of sDWI, up to the theoretical limits imposed by a mono-exponential signal model.Beyond b-values of ~ 1000 s/mm 2 , kurtosis effects stemming from the multi-compartmental nature of tissue start to become apparent in the diffusion signal decay, which cannot be accurately modeled using a simple mono-exponential function [33].Specifically, non-gaussian diffusion of restricted intracellular water results in more signal at high b-values than would be expected for purely gaussian diffusion [34], and therefore a systematic underestimation of high b-value signal by sDWI (Fig. 1).In this study, we showed that the median signal within the prostate was 55-72% lower on sDWI than aDWI.Ideally, b-values greater than 1000 s/mm 2 should be acquired to ensure that such non-gaussian effects are appropriately measured, enabling the assessment of restricted diffusion and derived biomarkers like RSIrs.There is evidence that measurements of restricted intracellular diffusion can help to differentiate between tumors of different histological patterns, like cribriform and non-cribriform, which differ in the degree of intracellular vs extracellular water [35][36][37].The clinical benefits of such granular assessment may offset the increase in scan time necessary to acquire data at higher b-values.Indeed, this study demonstrated a clear improvement in the detection of biopsy-proven csPCa with RSIrs compared to sDWI or aDWI.
Within the prostate and the prostate plus 5 mm margin, lesion conspicuity was reasonably preserved with sDWI.Both sDWI and aDWI had a similar quantitative performance in detecting csPCa with an AUC ranging between 0.56-0.65.However, sDWI introduced larger errors in the surrounding pelvic tissue even in a reduced FOV acquisition.Because the mono-exponential fitting was performed using a linear fit to log-transformed data for computational efficiency, areas with low SNR were susceptible to significant fitting errors.The pelvic region surrounding the prostate has many such regions with low SNR, including bone and connective tissue with inherently low SNR at the relatively long TEs used in this study, as well as fatty tissue with suppressed signal from the water-selective excitation pulse used during image acquisition.
Severe artifacts were observed on sDWI in these low-SNR regions of the pelvis surrounding the prostate, in particular for sDWI 500 , which makes the detection of metastasis outside of the prostate region difficult (Fig. 5).There are many ways to potentially improve the calculation of synthesized images including bi-exponential or multi-exponential modeling [23].For example, RSIrs is based on a multi-exponential model and may synthesize images without introducing artifacts.Image artifacts may be explained by poor signal quality, magnitude smaller than one in a subset of voxels, or noise/distortion correction post image acquisition leading to voxels with extremely low signal intensity.In particular, monoexponential models fail to correctly represent voxels with low signal intensity due to exponential fitting.Smoother images could be created by censoring those voxels by interpolating from surrounding voxels, smoothing low b-value images prior to calculation, or by thresholding low intensity voxels.For quantitative imaging, the details of such decisions would need to be clearly described and accounted for, and potentially could lead to more false positive/negative detections.Such enhanced images would not represent the measured truth and would include some unreliable voxels, which must be taken into consideration when interpreting the images.A "nicer looking" image does not necessarily mean that the image quality or reliability is better.
Prior studies reported sDWI to have higher subjective quality and tumor conspicuity [5,27].This may reflect the particular imaging sequences and, platforms used, or the particular image enhancement effects.It is also important to note that these prior results were mostly subjective judgements and not quantitative assessments of imaging quality.In the presented study we have proven that tumor conspicuity is quantitatively greater with aDWI (CNR = 0.95) in comparison to sDWI   One limiting factor of our study was that we only considered a retrospective dataset from a single scanner and a single institution.Also, the slice thickness of 6 mm for DWI data analyzed here is larger than the 4 mm recommended in PI-RADS v2.1.The larger through-plane voxel size increases signal-to-noise ratio but it is possible some very small tumors could have gone undetected, though any such inaccuracies would apply to aDWI, sDWI, and RSIrs alike.Further, only the conventional mono-exponential model was tested in the presented study, as this is the conventional method most cited and used for synthesis of high b-value DWI.A precise comparison of all possible methods for synthesizing DWI is beyond the scope of this manuscript.RSIrs is one quantitative biomarker based on a multicompartment model.The acquisition protocol in the datasets here was not optimized for models like hybrid multi-dimensional MRI or VERDICT.

Conclusions
Within the prostate, sDWI is a systematically inaccurate representation of aDWI, but the techniques are quantitatively comparable in terms of detecting csPCa with an AUC range between 0.56-0.65.In the surrounding pelvic tissue, high signal intensity artifacts are introduced with sDWI.These artifacts decrease CNR and thus affect the csPCa detection sensitivity in surrounding tissues and might mask potential metastases within the pelvis.RSIrs is superior to either sDWI or aDWI for quantitative csPCa detection.Despite the quantitative inaccuracies, sDWI may still be adequate for current subject clinical interpretation within the prostate.

Fig. 1
Fig. 1 Comparison of acquired images to those synthesized with mono-exponential models are presented for one representative patient.Mean values within the prostate using either b-values up to 500 s/mm 2 (sDWI 500 ) or b-values up to 1000 (sDWI 1000 ) are compared to acquired DWI (aDWI).are compared.Fig. A-F show the different diffusion images for one patient.A-D presents the acquired images for b = 0 s/mm 2 (A), b = 500 s/mm 2 (B), b = 1000 s/mm 2 (C) and b = 2000s/mm 2 (D).E and F show the synthesized b = 2000s/mm 2 images.E shows sDWI 500 and F sDWI 1000 .aDWI = acquired diffusion-weighted image for b = 2000s/mm. 2

Fig. 2
Fig. 2 Representative images from three patients (corresponding to three columns).A Acquired diffusion-weighted image (aDWI) for b = 2000s/ mm 2 , (B) synthesized DWI using acquired b-values up to b = 1000 s/mm 2 (sDWI 1000 ), (C) synthesized DWI using acquired values up to b = 500 s/mm 2 (sDWI 500 ), and (D) restriction spectrum imaging restriction score (RSIrs).The radiologist-defined cancer lesion for each patient is indicated in blue.All presented patients had a PI-RADS score of 5.The same window level was chosen for all presented images

Fig. 3
Fig. 3 Violin plots summarizing the signal intensity across 151 patients for (A) 50 th percentile, (B) 95 th percentile, (C) and 98 th percentiles of various DWI metrics calculated for each patient.The percentiles are estimated over different regions of interest: the prostate; the prostate with varying margin (5 mm, 30 mm, or 70 mm); and the whole field of view.aDWI = acquired diffusion-weighted image with b = 2000s/mm 2 ; sDWI = synthesized DWI for b = 2000s/mm 2 using either acquired b-values up to 1000 s/mm 2 (sDWI 1000 ) or up to 500 s/mm. 2 (sDWI 500 )

Table 1
Patient Characteristics range Q 1 -Q 3 Range between lower first quartile to upper third quartile, csPCa Clinically significant prostate cancer.MRI Magnet resonance imaging, PSA Prostate-specific antigen